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Abstract 

Fully reconstructed ti — > W + bW~b — ► Ivqq'bb events are used to deter- 
mine the fractions of right-handed (/+) and longitudinally polarized (/o) W 
bosons produced in top-quark decays. The helicity fractions are sensitive to 
the couplings and the Dirac structure of the Wib vertex. This paper reports 
measurements of the W-boson helicity fractions from two different methods 
using data corresponding to an integrated luminosity of 1.9 fb -1 of pp colli- 
sions at a center-of-mass energy of 1.96 TeV collected by the CDF II detec- 
tor operating at the Fcrmilab Tevatron. Combining the results from the two 
methods, we find /o = 0.62 ± 0.10 (stat) ± 0.05 (syst) under the assumption 
that /+ = 0, and /+ = -0.04 ± 0.04 (stat) ± 0.03 (syst) with f fixed to 
the theoretically expected value of 0.70. Model-independent fits are also per- 
formed and simultaneously determine /o = 0.66 ± 0.16 (stat) ± 0.05 (syst) and 
/ + = —0.03 ± 0.06 (stat) ± 0.03 (syst). All these results are consistent with 
standard model expectations. 
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1. Introduction 

Charged current weak interactions proceed via the exchange of a W ± boson 
and are theoretically described by a vertex factor that has a pure vector minus 
axial- vector (V — A) structure [1]. While weak interactions have been tested 
with high precision at low momentum transfers, e.g. in radioactive /3-decay, the 
vertex structure may be altered in interactions at high momentum transfers due 
to new physics contributions. Among the known fundamental particles, the top 
quark stands out as the heaviest, with a mass of m t = 172.4 ± 1.2 GeV/c 2 [2], 
and thereby gives access to high momentum scales. It has been suggested that 
the top quark may have non-universal gauge couplings as a result of dynamical 
breaking of the electroweak symmetry [3]. 

Given our present knowledge of the Cabibbo-Kobayashi-Maskawa quark- 
mixing matrix [4], the top quark decays with a branching ratio close to 100% in 
the mode t — ► bW + . The Dirac structure of the Wtb vertex can be generalized 
by the interaction Lagrangian 



r - in. 

V2 



W-lrf{fiP- + fiP+)t d v W-bo» v {tiP- + f*P+)t 



+ h.c. , 
(1) 

where P± = |(1 ± 7 5 ) and ia^ v = — §[7'*,7 I/ ] [5]. In general the interaction 
of fermions and gauge bosons can be expressed by six form factors. Assuming 
the W boson to be on-shell, the number of form factors is reduced to four. 
These four form factors f^'^ can assume complex values in general, but take 
values of /f = 1 and fi = f£ = = in standard electroweak theory, 
such that the production of right-handed W bosons from top-quark decay is 
suppressed. A general strategy to experimentally determine all four form factors 
in Eq. 1 involves the measurement of the VF-boson helicity fractions and the 
measurement of the single top-quark production cross-section in the ^-channel 
and in the s-channel [6] . 

The production of longitudinally polarized W bosons is enhanced due to the 
large Yukawa coupling of the top quark to the Higgs field responsible for elec- 
troweak symmetry breaking (EWSB). The fraction of right-handed W bosons, 
/+, is predicted to be very small O(10~ 4 ) [7], which is well below the sensitivity 
of the measurements reported here. The partial decay widths into the different 
W-boson helicity states explicitly depend on the form factors. Assuming the 
standard electroweak theory values for the form factors the fraction of longitu- 
dinally polarized W bosons is given by f = r{Wo)+ rlw-)+r(w+) ~ 2 m %\ra\ I 5 1 
at leading order in perturbation theory, where Wo and W± indicate longitudi- 
nally and transversely polarized W bosons, respectively. For m t as given above 
and a VF-boson mass of m w = 80.403 ± 0.029 GcV/c 2 [4] the theory predicts 
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fo = 0.697 ± 0.002. Next-to-leading-order corrections decrease the total decay 
width, as well as r(T / ^o), by about 10% [8], while f is only changed by about 
1% [9]. A significant deviation of fo or / + from the predictions exceeding the 
1% level would be a clear indication of new physics. 

This article reports the results of two analyses using the same dataset and 
their combination. Both analyses use the observable cos 9* , which is the cosine 
of the decay angle of the charged lepton in the VF-boson decay frame measured 
with respect to the top-quark direction. This has the following distribution: 

w(fl*) = /o-wo(n + /+-w+(e*) + (l-/o-/+)-w-(e*) with (2) 
= 7(1- cos 2 9*),uj+(9*) = !(i + cosr)V_(n = |(l-cosr) 2 . (3) 

4 8 8 

The parameters fo and /+ are the VF-boson helicity fractions to be determined. 

The two analyses estimate cos 9* for each event by reconstructing the full ti 
kinematics. These methods of reconstructing the four-vectors of the top-quark 
and antitop-quark as well as their decay products [10, 11, 12] possess a broad 
applicability and offer the possibility to measure a full set of top-quark proper- 
ties, such as the top-quark mass and the forward-backward charge asymmetry in 
ti production [13]. Experimental acceptances and resolutions introduce distor- 
tions of the cos 9* distribution which must also be taken into account. The two 
analyses employ alternative methods for reconstructing the ti kinematics, for 
correcting the experimental effects, and for determining the polarization frac- 
tions from the resulting cos 9* distributions in the observed events. They have 
similar sensitivities and are combined, taking into account correlations, to yield 
the most precise estimates of fo and / + . Both analyses subject the observed 
data to fits in three different scenarios: 

1. Measure fo under the assumption that /+ = 0. This corresponds to a 
model in which the form factors f^ and f% are zero, meaning there are 
no right-handed bottom-quark couplings present. 

2. Measure /+ under the assumption that fo = 0.7, which is sensitive to 
models with f% = f^ = 0, i.e. the presence of an additional V + A 
current in top-quark decay, but no additional magnetic couplings. Using 
the relation /+//- = {fi/fi) 2 one can translate the measured helicity 
fractions into the ratio of form factors. 

3. Measure fo and /+ simultaneously in a two-parameter fit, which is model- 
independent. 

Model-dependent measurements of fo and / + using smaller datasets have 
been previously reported by the CDF [14] and D0 [15] collaborations. Most 
recently the D0 collaboration has reported a model-independent result using 
1 fb _1 [15] of Tevatron data. The measurements reported here use twice as much 
data and improved analysis techniques and yield the most precise determinations 
of the M^-boson helicity fractions in top-quark decays. 
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2. Selection of tt candidate events 



The data used for the analyses reported here are collected by the CDF II de- 
tector [16]. We select events of the type tt — > W + bW~b — > tvqq'bb, which yield 
an experimental signature of one high energy charged lepton, missing trans- 
verse energy due to the undetected neutrino, and at least four jets, two of which 
are 6-quark jets. Exactly one isolated electron candidate with transverse en- 
ergy [17] Et > 20 GeV and pseudorapidity [17] < 1.1 is required, or exactly 
one isolated muon candidate with transverse momentum [17] Pt > 20 GeV/c 
and | r\ \ < 1.0. An electron or muon candidate is considered isolated if the Et 
not assigned to the lepton in a cone of R = \J (Arj) 2 + (A</>) 2 = 0.4, centered 
around the lepton, is less than 10% of the lepton Et or Pt, respectively. Jets 
are reconstructed by summing calorimeter energy in a cone of radius R = 0.4. 
The energy of the jets is corrected for differences as a function of rj, time, and 
additional energy depositions due to multiple interactions occurring in the same 
event [18]. An additional correction leads from calorimeter based jets to jets 
at the particle level. Candidate jets must have corrected Et > 20 GeV and 
detector |?7| < 2. Events are required to have at least four jets. The corrected 
missing transverse energy [19] $t accounts for the energy corrections made for 
all jets with corrected Et > 12 GeV and < 2.4 and for muons and is re- 
quired to be greater than 20 GeV. At least one jet in the event has to contain a 
secondary vertex identified using the algorithm described in [20] and consistent 
with having originated from a 6-hadron decay. Additional requirements further 
reduce the background contribution as follows. Electron events are rejected if 
the electrons originate from a conversion of a photon. Cosmic ray muon events 
are rejected as well. To remove Z bosons, events in which the charged lepton 
can be paired with any more loosely defined jet or lepton to form an invariant 
mass consistent with the Z peak, 76 — 106 GcV/c 2 , are excluded. With these 
selection criteria, we select 484 tt candidates in a sample corresponding to a 
total integrated luminosity of 1.9 fb _1 . 

Kinematic resolutions and selection and reconstruction efficiencies for tt 
events are determined utilizing pythia [21] and herwig [22] event generators 
where the top-quark mass is set to 175 GeV/c 2 . Samples of events generated 
with pythia, ALPGEN [23], and madevent [24], interfaced to pythia parton 
showering are used to determine certain background rates and to estimate the 
cos6>* distribution for background events. In order to develop and validate the 
methods presented, madevent and a custom version of HERWIG are used to 
generate samples of simulated events with controllable VF-boson helicity frac- 
tions. All generated events are passed through the CDF detector simulation [25] 
and then reconstructed in the same way as the observed events. 

3. Background Estimation 

The selected tt sample is estimated to be contaminated with about 87 events 
coming from background processes. These non-tt processes originate mainly 
from PF+jets events with a falsely reconstructed secondary vertex (Mistags), 
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from VF+jets events in which the jets are real b- and c-quark jets (VF+heavy 
flavor), and multi-jet processes that contain no real W boson (non-W). These 
backgrounds are estimated using a combination of data and Monte Carlo meth- 
ods as described in detail in [20] . Additional sources of background arise from 
electroweak processes like diboson production (WW, WZ, ZZ), the production 
of single top-quarks, and Z bosons. These backgrounds are predicted based 
on their theoretical cross sections and acceptances and efficiencies, which are 
derived from simulated events. Table 1 shows the background estimation and 
the observed number of events after all selection criteria. 



Background source 


N(>4jet) 


Fy+hcavy flavor 


37 ± 10 


Mistags 


20 ±5 


non-W 


18 ± 16 


Electroweak 


12 ± 1 


Total background 


87 ±23 


Observed events 


484 



Table 1: Expected number of background events and the number of observed events in 
a 1.9 fb _1 data sample using the selection criteria described in the text. 



4. Extraction of the W^-boson helicity fractions 

In order to measure the W^-boson helicity fractions we follow two approaches. 
Both analyses use cos 8* as the sensitive observable, estimated on an event-by- 
event basis by fully reconstructing the tt kinematics. The cos 9* distribution can 
be decomposed into three separate components according to the three different 
W-boson helicity states. The first analysis is based on the methods developed 
to precisely measure the top-quark mass [10] and uses the fact that the three 
helicity components have distinguishable shapes. In this technique we find the 
expected distributions ("templates") of the helicity components, containing res- 
olution effects, and superpose those. The helicity fractions are then given by 
normalizations from an unbinncd likelihood fit and the results are corrected 
for acceptance effects afterwards [26]. We refer to this analysis as the "tem- 
plate analysis" in the following. The second analysis, called the "convolution 
analysis", is based on the method described in [11, 12, 14]. Starting from the 
theoretically predicted number of events in each bin of the particle level cos 9* 
distribution we convolute acceptance and resolution effects with these predic- 
tions to derive the expected number of events in each bin of the reconstructed 
cos#* distribution. In this method, / and / + are then determined from a 
binned likelihood fit. 

The event selection and reconstruction of the two techniques employ different 
choices in the design of background suppression, jet flavor identification, and 
parton assignment. The agreement between the two methods shows that these 
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design choices do not bias the final result. While the convolution analysis uses 
the standard event selection described in sec. 2, the template analysis chooses 
to place an additional cut on the scalar sum of all transverse energies of the 
event, and requires Ht > 250 GeV to further suppress multi-jet non-W 
background. This results in 53 ± 20 events estimated as background, and 
reduces the total number of selected events to 430. A combinatoric ambiguity 
arises in the reconstruction of the ti kinematics when choosing which of the 
reconstructed jets corresponds to which of the final state quarks in the ti — > 
tvqq'bb decay. The analyses each test all possible jet-quark assignments and then 
use alternative criteria to choose the "best" one for each event. The template 
analysis uses the technique described in [10]: jet energies float within expected 
resolutions, 6-tagged jets are assigned to b quarks, and the top-quark mass is 
left floating in the fit while the IT-boson masses are constrained to 80.4 GeV/c 2 . 
The algorithm described in [11, 12, 14] is used in the convolution analysis. The 
jet-quark assignment is selected using constraints on the IF-boson mass, the ti 
mass difference, the transverse energy in the reconstructed ti pair with respect 
to the total transverse energy in the event, and the &-jet probability of the 
jets. Neither analysis assumes a particular value for m t in the reconstruction; 
since /o has an explicit mj-dependence, doing so would introduce a bias in the 
measurement. Although the algorithms to reconstruct the kinematics of the ti 
pairs are different, the cos 6* resolution for each analysis is estimated to be the 
same (~ 0.35) from studies using generated ti events. 

In both analyses the VT-boson helicity fractions are determined from max- 
imum likelihood fits to the resulting cos#* distributions. The two analyses 
employ alternative methods to derive the fit inputs which will be discussed in 
more detail in the next paragraphs. In the fits, the helicity fractions /o and 
/ + are free parameters, the constraint /_ = (1 — /o — /+) is applied, and the 
background contribution is allowed to float but is Gaussian constrained using 
an RMS corresponding to the uncertainty on the estimate of the total number 
of background events. As already discussed in sec. 1, each analysis performs 
three different measurements. In two measurements we determine fo or / + and 
fix the other parameter to the value expected in case of a pure V — A structure 
of the Wtb vertex (/ = 0.7, /+ = 0.0). In the third measurement, f and / + 
are both treated as free parameters and are measured simultaneously. 

The template method utilizes samples of generated ti events in which the 
leptonically decaying W boson is forced to a specific polarization to get the 
normalized cos 9* probability distribution function V(cos9*) for each FT-boson 
polarization. These generated events satisfy all the selection criteria and are 
reconstructed in the same manner as the observed events. The P(cos9*) for a 
certain helicity mode is obtained by fitting the reconstructed cos 9* distribution 
obtained from the corresponding generated ti events and does not depend on 
the helicity fractions assumed for the hadronically decaying W boson. The 
background modeling is verified by comparing the distribution obtained from 
generated events to the distribution of observed events in which there is no 
secondary vertex tag and in those for which the decay length of the secondary 
vertex tag is negative, meaning that the reconstructed secondary vertex and 
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Figure 1: The V(cos9*) used in the template analysis, which are the reconstructed 
ui(9*) distributions for longitudinal, right- and left-handed W-boson helicities, as well 
as the u(8*) in the background model. The curves are normalized to the same area. 



the reconstructed jet itself are located in opposite hemispheres with respect to 
the primary vertex. These are background dominated samples. The V(cos9*) 
parameterizations are empirically chosen to provide a good description of the 
cos#* distributions and use a third degree polynomial times two exponential 
functions. The resulting V{cos9*) are compared in Fig. 1. Using alternative fit 
functions, negligibly affects the results. 

Since the kinematics of the VF-boson decay depend on its polarization, the 
kinematic cuts applied have different acceptances for the different polarizations 
and alter the observed composition of polarization states. The largest impact 
is due to the isolation requirement and the cut on the px of the charged lepton. 
Therefore a correction is applied to the obtained helicity fractions to account 
for these acceptance effects before presenting the results. 

In the convolution analysis the cos 9* distribution is reconstructed in six 
bins, corresponding to the resolution of the reconstruction of the it kinematics. 
The starting point for the extraction of the VF-boson helicity fractions in this 
method is the theoretically predicted number of signal events in each bin of 
the cosO* distribution, // lg (/n, /+), depending on /n and /+, which can be 
calculated using Eq. 2. Acceptance and resolution effects are then taken into 
account [14] by convoluting both effects with the theory prediction. This leads to 
the number of signal events expected to be observed in a certain bin accounting 
for all distorting effects: 



The migration matrix element S(i, k) gives the probability for an event which 
was generated in bin i to occur in bin k of the reconstructed cos 9* distribution. 
Since the acceptance depends on cos 9* , we weight the contribution of each bin 
with its event selection efficiency e^. The effects considered are independent 
of the W-boson helicity fractions and this is validated using several samples 
of generated events with different VF-boson polarizations. Thus, and S(i,k) 



/4 is '° bs (/o,/ + ) « 5>f S (jb,/+) . £j . S(i,k). 



(4) 
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source 


Sfo 




s.U 




Sfo 




Sf+ 




f + fixed 


fo fixed 


combined fit 


combined fit 




tempi. 


conv. 


tempi. 


conv. 


tempi. 


conv. 


tempi. 


conv. 


JES 


0.024 


0.045 


0.017 


0.025 


0.021 


0.016 


0.027 


0.032 


ISR 


0.002 


0.010 


0.003 


0.003 


0.010 


0.036 


0.007 


0.014 


FSR 


0.021 


0.025 


0.009 


0.011 


0.025 


0.045 


0.002 


0.016 


Bkg 


0.023 


0.032 


0.016 


0.019 


0.018 


0.028 


0.017 


0.032 


MC 


0.019 


0.012 


0.009 


0.005 


0.019 


0.015 


0.010 


0.002 


PDF 


0.005 


0.005 


0.005 


0.002 


0.005 


0.014 


0.002 


0.006 


Total 


0.044 


0.062 


0.027 


0.034 


0.043 


0.072 


0.034 


0.050 



Table 2: The sources of systematic uncertainties and their related estimates for the 
template analysis (tempi.) and the convolution analysis (conv.). The total systematic 
uncertainty is taken as the quadrature sum of the individual sources. 



can be estimated from a sample of events generated with the pythia event 
generator using the standard settings. The total number of events expected to 
be observed in a certain bin is then given by the sum of M^ ls ' obs (/o, /+) and the 
expected number of background events, which is independent of the IT-boson 
polarization and is derived from the background composition shown in Table 1. 
In a maximumdikclihood fit the expected number of events is compared bin by 
bin to the number of observed events to determine fo and / + . 

In order to compare our observations with theory, we subtract the back- 
ground estimate from the reconstructed cos#* distribution, correct for accep- 
tance and resolution, and normalize the distribution to the ti cross section of 
<j t t = 6.7 ± 0.9 pb [27, 28]. The correction is made by applying a bin-by-bin 
correction factor to the cos#* distribution. The correction factor is given by 
/4 ig (/o t ,/+) divided by /iJ ig,obB (/^ t , ff), where / fit and ff are the obtained 
results. 

The systematic uncertainties associated with the measurement of fo and /+ 
are summarized in Table 2. The systematic uncertainties were determined by 
constructing ensemble tests with signal and/or background templates, affected 
by the systematic under study, but fit using the default parameterizations and 
normalizations described above. We studied the influence of variations in the 
jet energy scale (JES) and of variations in initial and final state radiation (ISR, 
FSR). The latter was estimated by producing samples of simulated events for 
which the simulation was altered to produce either less or more gluon radiation 
compared to the standard setting [10]. Specifically, two parameters controlling 
the parton shower in the pythia program are varied: Aqcd and the scale factor 
K to the transverse momentum scale of the showering. The different settings are 
derived from studies of ISR in Drell-Yan events. We also studied the influence 
of the background modeling (Bkg), of different Monte Carlo event generators 
(MC), and of the parton distribution function (PDF). The resulting shifts in 
the mean fitted longitudinal and right-handed fraction are used to quantify 
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cos 8 

Figure 2: The observed cos 9* distribution (points) overlaid with the fit-curves for the 
three different fit-scenarios (as explained in sec. 1) for the template analysis. 



the systematic uncertainties. The positive and negative variations obtained are 
symmetrized by choosing the maximum deviation. The ensemble tests were all 
performed using pythia generated events with m t = 175 GeV/c 2 as signal with 
the W-boson helicity fractions fo = 0.70 and / + = 0.0, and the background 
model as described above. We have verified that these uncertainties do not 
depend on the actual value of f and / + by fitting samples of generated events 
with different VF-boson polarizations. 

The analyses presented in this paper use a top-quark mass of 175 GeV/c 2 . 
Since fo explicitly depends on the top-quark mass, the dependency of the mea- 
sured value of fo on the top quark mass is not treated as a systematic uncer- 
tainty. The measured value of /+ is only negligibly affected by variations in the 
assumed top-quark mass. 

5. Results and Combination of the Results 

The cos 9* distribution from the observed events is shown in Fig. 2 and 
Fig. 3 for both analyses together with the fits for fo and /+ and the model 
independent measurement. The results for the three different measurements 
together with the statistical and systematical uncertainties in both analyses are 
summarized in Table 3. In the template analysis the correlation between fo and 
/+ is determined to be -0.87 in the simultaneous fit, while for the convolution 
analysis the correlation is -0.89. 

We combine the single results accounting for correlations using the BLUE 
method [29]. The combined results can be found in Table 3. The statistical cor- 
relation between both analyses is estimated from ensemble tests using samples 
of generated events which account for the event overlap in the signal contri- 
bution. For the two model-dependent scenarios the correlation coefficients are 
found to be 0.66 and 0.65 when fitting for fo or /+, respectively. The correlation 
matrix for the model-independent scenario is given in Table 4. The resulting 
combination is weighted towards the template determination of /+ since its total 
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Figure 3: On the left-hand side the observed cos 6* distribution (points) is presented 
overlaid with the fits for /o and /+ for the convolution analysis. On the right-hand side 
the deconvoluted (using the fit result of the /+ measurement) distribution normalized 
to the tt cross-section is shown together with the theoretically predicted curves for 
purely left-handed, right-handed, and longitudinally polarized W bosons. 



uncertainty is significantly smaller than the total uncertainty from the convolu- 
tion method. Due to the strong anti-correlation between /o and /+ (see Table 4) 
the fo determination is affected correspondingly. The systematic uncertainties 
are taken to be completely correlated between the two methods. When combin- 
ing the model-independent results the systematic uncertainties for /o and /+ 
are taken to be 100% anti-correlated. The combined values of /o and / + have 
a correlation of —0.82. The combination improves the sensitivity by about 10% 
relative to the measurements of either method separately. 

In conclusion, we present two different analyses and their combination de- 
termining the VT-boson helicity fractions in top-quark decays, giving the world's 
most sensitive result for measuring these fractions so far. In addition to mea- 
suring fo and /+ separately, while fixing the other parameter to its expected 
value, we present a model-independent simultaneous measurement of the two 
fractions. All of these results are consistent with the values predicted within 
the electroweak theory of the Wtb vertex. 







template 


convolution 


combination 


X'/dof 


M/+ 


= 0.0) 


0.59 ±0.11 ±0.04 


0.66 ±0.10 ±0.06 


0.62 ±0.10 ±0.05 


0.7/1 


f+(fo 


= 0.7) 


-0.04 ±0.04 ±0.03 


0.01 ±0.05 ±0.03 


-0.04 ±0.04 ±0.03 


1.8/1 


fo 




0.65 ±0.19 ±0.04 


0.38 ±0.21 ±0.07 


0.66 ±0.16 ±0.05 


4.3/2 


u 




-0.03 ±0.07 ±0.03 


0.15 ±0.10 ±0.05 


-0.03 ±0.06 ±0.03 


4.3/2 



Table 3: Results of the template analysis, the convolution analysis, and the combined 
values. The results are given together with their statistical and systematical uncer- 
tainties. In addition the x 2 /dof of the combination is given. 
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template f 


convolution fo 


template / + 


convolution / + 


template /o 


1.00 


0.45 


-0.87 


-0.40 


convolution / 


0.45 


1.00 


-0.42 


-0.89 


template / + 


-0.87 


-0.42 


1.00 


0.48 


convolution /+ 


-0.40 


-0.89 


0.48 


1.00 



Table 4: Correlation matrix for combining the template and convolution analyses in 
the model-independent scenario. 
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